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Abstract 

Vilkovisky has claimed to have solved the black hole backreaction problem 
and finds that black holes lose only ten percent of their mass to Hawking radi- 
ation before evaporation ceases. We examine the implications of this scenario 
for cold dark matter, assuming that primordial black holes are created during 
the reheating period after inflation. The mass spectrum is expected to be domi- 
nated by 10-gram black holes. Nucleosynthesis constraints and the requirement 
that the earth presently exist do not come close to ruling out such black holes as 
dark matter candidates. They also evade the demand that the photon density 
produced by evaporating primordial black holes does not exceed the present 
cosmic radiation background by a factor of about one thousand. 

PACS: 4.70.-S, 4.70.Dy,95.35.+d,98.80.-k.,98.80.Ft 

Keywords: Cosmology, primordial black holes, dark matter, Hawking radiation, 
evaporation. 

1 Quasi-evaporating black holes 

In a recent series of papers, Vilkovisky [1, 2, 3, 4] has claimed to have completely 
solved the black hole backreaction problem, which is the major unresolved issue re- 
garding black hole physics since Hawking's discovery of black hole radiation thirty 
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years ago [5]. Hawking's original calculation was semiclassical and ignored the effect 
of the black hole's evaporation on the spacetime metric itself; Hawking argued that 
any such backreaction should be negligible until the black hole approaches the Planck 
mass. Vilkovisky maintains that if backreaction is properly taken into account, it ex- 
erts a negative feedback on the evaporation process, which ceases after the hole has 
lost only ten percent of its mass to thermal radiation. If the scenario is correct, then 
contrary to the accepted picture, primordial black holes of mass less than 10 15 g will 
not have entirely disappeared by the present epoch of the universe and could conceiv- 
ably provide the cold dark matter sought by cosmologists. Vilkovisky's calculations 
are also semiclassical and do not assume any exotic physics. Nevertheless, they have 
yet to receive full scrutiny by the physics community and it is premature to claim 
they are correct. In this paper we merely propose to take his scenario seriously and 
examine the consequences. We find that, indeed, black holes with mass less than 10 15 
grams can provide the Q ~ .3 dark matter while safely evading basic observational 
constraints. Such a conclusion may seem fairly obvious, since the black holes are 
evaporating only ten percent of their mass. On the other hand, their evaporation 
profile does not at all resemble the usual one and it is worth verifying that this does 
not seriously alter the expected conclusions. Furthermore, other phenomena, such as 
vortons[6], behave similarly at least insofar as that they are stabilized against contrac- 
tion by the current they generate. We do also feel that it is worthwhile to stimulate 
discussion of Vilkovisky's work, which if correct, may have provided the answer to 
two outstanding problems. 

Carr et al. have provided several reviews of the standard picture of primordial 
black hole (pbh) formation processes and observational constraints on their mass den- 
sity; see for example [7, 8, 9] and references therein. Although the density of black 
holes with M > 10 15 grams is constrained by a variety of techniques, including mi- 
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crolensing and 7-ray bursts, all constraints on holes with M < 10 15 grams are imposed 
by assuming that they interfere with established cosmological processes through their 
evaporation, which in the Hawking picture ends in an explosive phase. In the Vilko- 
visky picture, by contrast, the temperature profile of the emitted radiation differs 
sharply from the usual one, being nearly flat until it drops to zero, and evaporation 
never reaches an explosive phase. 

The Hawking temperature is the usual one, 



where M p i is the Planck mass, k is the surface gravity and u = t — r is the retarded 
time. (We work in units with c = G = ksoitzmann — 1-) However, in Vilkovisky's 
scenario the surface gravity is not k — His considerations of the radiation at 
infinity lead him to describe the metric in terms of two "Bondi charges." The first 
is Ai, the Bondi mass, which is the mass remaining in a compact region after the 
escape of radiation (the ADM mass minus the mass of the radiation). The second is 
a "charge" Q, which the black hole manufactures from the vacuum to protect itself 
from quantum evaporation. We emphasize, however, that these two "charges" are 
not merely assumed to exist but are a necessary consequence of the radiation metric. 
They are in fact, the two leading coefficients of the metric expansion at infinity and are 
the only two that need to be specified. Because in terms of these two coefficients, the 
asymptotic metric assumes the Reissner- Nordstrom form, the surface gravity becomes 
that for a Reissner-Nordstrom black hole: 

(-^ 2 -Q 2 ) 1/2 (12) 
[M + (M 2 - Q 2 )V2]2 ^ • > 

The evolution of M. and Q are governed by two coupled differential equations: 



dM M. 



2 

p 1 .Ji 



k (1.3) 



du 487r 

IT - £«■ (1-4) 
du 24vr v ' 



Quasi-evaporating black holes... 4 

We plot A4, Q and T in figure 1. As one sees, T is essentially flat but reaches a 
slight maximum of about 

T max = — , (1.5) 

before dropping abruptly to zero, at which point evaporation stops. Here M is the 
initial mass of the black hole. 
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Figure 1: M., Q 2 and T plotted in terms of the initial mass vs. normalized retarded time u for the 
Vilkovisky scenario. Notice that the temperature is not proportional to 1/M. but is reasonably flat 
until it drops to zero and that M. has dropped to about .9 of its initial value by the time evaporation 
ceases. 
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where M is in grams. 

The retarded time taken for evaporation is 

M 3 



Am « 96tt 



(1.6) 



(1.7) 



'U 



where u is evaluated at X + . We assume that asymptotically u becomes the cosmic 
time in a flat Friedmann-Lemaitre-Robertson- Walker universe, which in the matter- 
dominated case yields Au = At(l — 3 (to/ At) 2 / 3 ) for initial time t . We see that 
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for At >> to we may take Am ps At. Thus the cosmic time over which the hole is 
evaporating 

M 3 

T evap ~ 96^—^ « 1.6 x 1(T 27 M 3 (1.9) 

for M in grams. 1 

Equating r e „ ap with the age of the universe t u ps 4.5 x 10 17 s, gives a value, 

M u ps 6.6 x 10 14 grams. (1.10) 

This is essentially the canonical 10 15 gram black hole but in Vilkovisky's scenario it 
is not the mass of a hole that is evaporating during the current epoch; it is merely 
the mass of a hole whose evaporation is currently ceasing, by which time it will have 
lost only about ten percent of its mass, as can be seen in figure 1. Thus, M u is 
the value above which black holes effectively do nothing whatsoever by the current 
age of the universe. We expect that such a change should have drastic astrophysical 
consequences. 

2 PBH spectrum and total mass constraints 

We would like to check whether, given a reasonable pbh mass spectrum, there are 
any obvious reasons that pbhs of M < 10 17 grams (where femtolensing constraints 
kick in [7]) could not comprise the dark matter density VL dm ps .3. There are both 
observational and theoretical reasons to believe that the density fluctuation spectrum 
in the early universe should be scale-invariant, and in fact only in this case can pbhs 
have an extended mass spectrum. Nevertheless, one expects any black holes formed 
before inflation to be diluted to negligible density by the exponential expansion and 

*As in Hawking's original work, Vilkovisky's calculation was carried out for a single scalar field. 
If N species of particles are emitted with equal probability, the evaporation time will decrease by 
a factor 1/iV; however Page's results[10] show that emission is highly suppressed for particles with 
nonzero spin and thus only spin-zero need to be considered. Unfortunately, we do not know the 
number of spin-0 species the energies of interest in this paper. If N < 10, the results to be discussed 
go through unchanged. If N > 10, the limits to be discussed will necessarily be relaxed even further, 
because the integration time over pbh evaporation will be drastically reduced. 
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that therefore any relict pbhs would have been formed after inflation, presumably 
during the reheating phase. 

In this case one faces a serious difficulty: Because reheating takes place as inflation 
ends, any wave mode perturbing the background at this time will have not entered 
the Hubble radius ("horizon") at the onset of inflation to be frozen in as it exited 
later on. Consequently there is no clear reason for the perturbation spectrum to be 
scale invariant. Carr[9, 7] gives the following expression for the probability that a 
spherical region of mass M forms a black hole: 



where p = 7p is the equation of state and M* is a fiducial mass. The prefactor 
is actually the ratio S/5 max , where 5 = dM/M is assumed to obey the power law 
S = £(M/M*)~ n and 5 max = (M/M ) _2//3 results in maximal black hole formation. 
In this formulation n = 2/3 represents the scale-invariant case and for n > 2/3 
significant pbh production can take place only when e ~ 1. From Eq. (2.1) one can 
show that for n = 2/3 the number density of pbhs today with mass between M and 
M + dM should be 



Here, p* is the radiation density at the epoch when black holes of mass M* were 
formed, Rq is the present-day scale factor and is the scale factor at the formation 
epoch, which we will assume to be the reheating period. The factor j3 = e exp ^— ^ j 
represents the fraction of regions of mass M that collapses into black holes. 
For n 7^ 2/3 one can show with somewhat more difficulty that 




(2.1) 




(2.2) 



dn(M) 




(2.3) 
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The n = 2/3 case is certainly the easiest spectrum to deal with and because it is a 
power-law, we expect it to be dominated by pbhs of M — M*. However, as mentioned, 
it is not obvious that pbhs forming as inflation ends should obey the scale-invariant 
spectrum. From the usual Fourier decomposition of density perturbations [11] one has 
5 = 5p/p oc k z l 2 bk with S k oc k a l 2 . Evaluated at the wavelength corresponding to the 
mass, A ~ k~ l ~ M 1 / 3 and (Sp/p) M = dM/M oc M~ {1 / 2+a l 6 \ Note that a = 1 gives 
5 oc M" 2//3 , the Harrison-Zel'dovich spectrum. From the Mukhanov equation[12], one 
expects that any modes that do not cross the Hubble radius during inflation should 
have a = 3. Equating exponents to the pbh- formation case, where 5 = (M / 'M*)~ n , 
gives n = a/6 + 1/2. For a = 3, we have n — \. 

Notice from Eq.(2.3), however, that for any n ^ 2/3, there is not only an exponen- 
tial cutoff but that dn(M) strictly vanishes when the last bracket equals zero. Since 
one expects e ~ 1 and 7 = 1/3 for a hard equation of state, the n = 1 case gives 
M max ~ 30M*. The exponential cutoff ensures that the number density will be negli- 
gible at an even smaller value of M max . From these considerations, it is reasonable to 
limit ourselves to the scale-invariant spectrum (for simplicity and ease of comparison) 
and to delta-function spectrum, in which all the pbh mass is concentrated in black 
holes of mass M*. As we will see, the results do not differ significantly. 

For the scale-invariant case with a radiation equation of state, Eq. (2.2) yields a 
total present mass density [9, 7] 

ml 12 \ M-^dM. (2.4) 
Since p rad oc i?~ 4 and p p bh oc R~ 3 , Eq. (2.4) can be rewritten as 



/•Mmax 

%bhpcrit = VradPcrit(l + Z*)f3Ml' 2 / M~ 3/2 dM, 



— 2f2 rac ;p cr .j t (l + z*)/3 



M. ^ 1/2 



(2.5) 



for formation redshift z*. Note that placing constraints on the mass density of pbhs 
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is equivalent to constraining (3. 

In the standard scenario, one typically cuts the spectrum off at the lower limit 
M* PS 10 15 grams because smaller holes would have evaporated by today. Taking the 
upper limit to be oo gives the total mass density above the cutoff. In Vilkovisky's 
scenario, as already explained, smaller pbhs are still present. If we take M* to be the 
minimum mass of pbhs formed after inflation, then we expect it to be roughly the 
horizon mass at the time of formation [7], or in a radiation-dominated universe 

M, = M pl {^j (2.6) 

where T rh is the reheating temperature. Since T rh rs 10 16 GeV and T pl ps 10 19 Gev, 
we have M* ps l0 6 M pt ps 10 g. Then, because M* <<< M max , Eq. (2.5) becomes 
simply 

tt pbh & VL rad (l + z,)(3, (2.7) 

independent of the upper limit M max . For fl pbh < .3, VL rad ~ 10~ 4 and ~ 10 29 , the 
above yields (3 ~ 10 -26 . 

To obtain a delta-function spectrum, we replace the M~ 2 dM in Eq. (2.3) by 
M~ 1 5(M — M*)dM, where now 5 is a Dirac delta function, and integrate over M, 
giving 

n = ^{^j M;\ (2.8) 

In this case /3 is a different functional form of e than the scale-invariant spectrum but 
the physical interpretation remains the same, the fraction of the volume of space that 
has collapsed into black holes. With the approximations we have made, Eq. (2.7) 
remains unchanged for the delta-function spectrum, and the limit on (3 the same. 
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3 Entropy and annihilation constraints 

The above limit was set merely by requiring the total mass density in primordial black 
holes not exceed the known dark matter density and it is not immediately obvious that 
this value for f3 is free of observational contradictions. Specifically, it does not take 
into account pbh evaporation, which can be expected to affect the cosmic background 
photon density. The strongest constraints [8] in the standard model assume that 
10 15 g pbhs are currently evaporating in the 7-region. Requiring that density of 7's 
from pbhs be less than that of the background, pushes (5 down to about 10" 28 . We 
now perform a similar analysis for the Vilkovisky scenario, although here there is a 
significant difference: Since black holes < 10 15 g continue to evaporate to the present 
at an approximately constant rate, for the scale-invariant spectrum we must integrate 
the photon emission over the entire history of the universe. (From Eq. (1.6) it is clear 
that the temperature of small black holes can be high enough to evaporate massive 
particles. In the present analysis, we assume that only massless particles, which we 
refer to as "photons," are emitted; this will only tighten the constraints. See also 
footnote 1.) 

For the scale-invariant case, if .1M worth of photons is evaporated per hole, then 
Eq. (2.2) allows us to write for the photon energy density due to pbhs at emission 



Here the subscript e denotes "at emission" and the temperature T oc R" 1 refers to 
the background temperature. This equation assumes that all pbhs are continuously 
evaporating, but after a time r evap given by Eq. (1.9), a black hole shuts off and no 
longer contributes to p 1 . Inverting Eq. (1.9) therefore gives in a lower-limit cutoff for 
M as a function of t. If is the smallest allowable cutoff and we take T oc t~ 2 / 3 
over most of the history of the universe, then M cut = M*{T T */T e ) l l 2 , where T T * is 
the temperature at the time when a hole of mass M* ceases to evaporate. Given a 




(3.1) 
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formation time t* for pbhs of mass M*, we have £*/t* = M^J (967rM 2 ). Assuming 
that formation takes place in a radiation-dominated phase, then T T */T* = (t^/V*) 1 / 2 
and thus T T * = (M p^ /^/967^M*)T*. Consequently, 

m -(^) 1/2m - 1/2 (I) 1/2 - 10 " 3m - 1/2 (I) 1/2 ' 

where the last expression is for M p i in grams. 

Now integrating Eq. (3.1) between M cut and M and retaining only the lower limit 
gives 

frp \ 13/4 

p 7e « .2 x 10 3 / 2 /3p,My 4 I ^\ . (3.2) 
Redshifted to the present background temperature T , p 7e becomes p 70 = p 7e (T /T e ) 4 



or 



p 70 « .2 x 10 3 / 2 /3p*My 4 (-g^) T" 3 / 4 . (3.3) 



This can be regarded as the integral of dp l0 evaluated at T e and T , or 

fm \ 13/4 /T x 3/4 

p 70 « .2 x 10 3 / 2 /3p«My 4 ^ j = .2 x 10 3 / 2 /3p My 4 J , (3.4) 

where we have once more discarded the upper limit and have set p* = po(T^,/T ) 4 , 
where po is the current radiation density. 

Since TJT ~ 10 29 and M* « 10 g, then p 70 ~ .3 x 10 23 - 5 /3p . If we require that 
p 7 o < po, then (5 < 3 x 1CT 23 - 5 . Because this is significantly larger than the /3 ~ 1CT 26 
found by requiring f2 p &/i < .3, 10-g pbhs are not ruled out as dark matter candidates 
at the level of this calculation for the scale-invariant spectrum. 

For the delta- function spectrum, we take p 7e = p p bh x (number of photon's per hole) 
.lMppbh/T^h, where T 7 /j is the energy of a 7 emitted by the black hole. Then from 
Eq. (2.8) 

TA 3 M, 



p 7e = .1(3 p* ( — J — , (3.5) 



Quasi-evaporating black holes... 



11 



or in terms of po an d redshifted to the present as above 




(3.6) 



We approximate T e by T T *, the background temperature at evaporation; as calculated 
above T r * ~ 10" 10 T p/ . From Eq. (1.7), T lh ~ 10^ 8 T p/ » T T * but photons from 
10-gram black holes will thermalize on a timescale far shorter than the expansion 
timescale and merely become an indistinguishable part of the background radiation 
density. Thus instead of T lh we use T T * in Eq. (3.6) and find that f3 ~ 10~ 22 , which 
is reasonably consistent with the scale-invariant calculation. 

We note that the constraint p 70 < po is the weakest that one can apply, but 
seems to be the appropriate one for these models. As just mentioned, any photon 
emitted before decoupling will thermalize. From Eq. (1.9) only photons from pbhs 
of M > 10 13 g will be emitted after decoupling and not thermalize. Eq. (1.6) shows 
that the energy of such photons is ~ 100 MeV, or about 100 KeV when redshifted to 
the present. The x-ray background in this part of the spectrum is indeed only ~ 10~ 5 
of the CMB background, but with the scale-invariant mass spectrum, holes in the 
10 13 -g range form a fraction only ~ 10~ 30 of the total and therefore can make only 
a negligible contribution to the x-ray background. For the delta-function spectrum 
such holes do not even exist. Nevertheless, any constraint that is a factor of ~ 10 3 
tighter than the one we have imposed will rule out these models. 

In searching for further constraints we should consider the following: Because pbhs 
of M ~ 10 g dominate either mass spectrum, for Vt p bh ~ 1 we must have n ~ 10~ 30 
cm -3 . Assuming the same density enhancement within the galaxy above the average 
density of the universe implies n ~ 10~ 25 . In that case collisions with the earth would 
be expected on a timescale r ~ (crnv ) _1 ~ 10~ 3 (c/v) s. For v ~ 200 km/s, consistent 
with galactic rotation curves, we have r ~ 1 s. Thus, over its lifetime, the earth 
has experienced ~ 10 17 black hole collisions. This is a somewhat alarming prospect; 
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one is therefore prompted to investigate the effects of accretion resulting from such 
interactions. 

In the simplest model[13, 14] (nonrelativistic, good in any number of dimensions), 
a sub-atomic-size black hole swallows anything that comes within a collision cross 
section: 



where r n ~ 10~ 13 cm is the radius of the nucleon, p is the average density of the 
medium and v is the black hole velocity. This accretion law is valid for r bh < r n and 
with v = constant yields trivially 



We may require that the earth presently exists. Since M e ~ 10 g, if each of the 
10 17 pbhs grew to 10 11 g, this would presumably be an indication that the earth was 
annihilated by accretion. The shortest accretion timescale r acc will be given by the 
maximum v. Taking v to be the impact velocity ~ 10 7 cm/s, and p e ~ 5 g/cm 3 , Eq. 
(3.8) however yields r acc ~ 10 29 s. 2 

More sophisticated analysis will not lower this number into the realm of measur- 
ability, regardless of spectrum. Hence, the continued existence of the earth puts no 
limits on the density of 10 g pbhs. Because r acc scales as M/a oc M 1//3 for an astro- 
nomical body of mass M, neither will observing the potential destruction of asteroids 
significantly improve the annihilation constraint. 

4 Nucleosynthesis Constraints 

Big Bang nucleosynthesis (BBN) has also been used to limit pbh density (see [15, 

16, 17, 18] and references therein). All BBN constraints on pbhs involve black holes 

2 The above accretion law can be written dM/dx ~ op or dM <~ m n dx/\, which states that the 
black hole picks up a nucleon's worth of mass every mean free path A, independent of velocity. Thus 
for a black hole to accrete Mf <~ 10 11 g requires 10 11 nucleons, or about 10 35 cm. 



dM 




(3.7) 



dt 



Mf ~ irr^pvt. 



(3.8) 



Quasi-evaporating black holes. .. 13 

of mass ~ 10 9-10 g, which are evaporating during helium production at 1 s< t < 
100 s. Various processes can be considered, which either spallate the helium as it 
is being produced and overproduce deuterium, or destroy the deuterium itself (by 
photodisassociation). For such scenarios the delta-function spectrum of 10-gram black 
holes is clearly inapplicable, since all the black holes will have ceased evaporation long 
before BBN, and because we have assumed that the emitted photons are thermalized, 
by BBN they can no longer be energetic enough to have any particular effect. 

For the scale-invariant mass spectrum, one should put in a lower mass cutoff, as we 
did in §3, but one appropriate for a radiation-dominated universe, where T oc t -1 / 2 : 

M cut « 10" 3 My 2 

However this cutoff implies that the only black holes continuing to evaporate at BBN 
are those with mass > 10 10 g, as expected. Thus we are effectively constraining black 
holes of same mass range as the earlier studies and may merely cite their results, all 
of which give /3 < 10~ 21-22 . These figures should be relaxed by an order of magnitude 
to reflect the fact that in the Vilkovisky scenario only AM is released as radiation. If 
we regard 10 9 -g black holes as the tail-end of the power-law spectrum, then a single 
value of P < 10~ 20 applies across the entire mass range. This is a very weak limit 
compared to that obtained by constraining the present-day density of pbh photons 
and does nothing to rule out 10-g pbhs as dark matter candidates. 

Thus, to conclude, the standard methods employed to limit pbh density place no 
limits on the black holes in the Vilkovisky picture. To some extent this was expected, 
since quasi-evaporating black holes lose only ten percent of their masses to radiation. 
In any case, although one might have theoretical prejudices against primordial black 
holes, such prejudices are not sufficient to rule them out as cold dark matter candi- 
dates. 
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